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Sputtering  of  Single  Crystals: 
Experimental  Evidence  of  the  Ejection  Process 

By  Nicholas  Winc^ad 
DepaftmcBt  of  Ckemittry,  Pena  State  Uaivenity 
152  Davey  Laboratory,  Uaivenity  Park,  PA  16802,  USA 


Synopsis 

The  objective  of  thie  review  k  to  munmatite  recent  experknental  evidence  of  the  eputteriiig  of 
atons  from  einsle  cryttal  targeta.  The  information  k  extracted  primarily  from  literature  pub- 
Ikhcd  within  the  laat  10  yean  and  k  meant  to  cam|diment  other  recent  tevkwa  on  lingk  cryetal 
sputtering.  The  emphaak  of  the  discuaaion  k  on  experiments  that  are  performed  on  wdl  diar- 
acterised  surfaces  obtained  under  low-dose  conditions.  Because  of  thk  restriction,  most  of  the 
experiments  utilised  laser  poet  ionisation  as  the  detection  scheme.  The  dkcnssion  of  experimen¬ 
tal  results  k  focused  mainly  toward  comparison  with  results  of  mdecular  dynamics  computer 
aintilations.  The  review  includea  a  discussion  of  experimental  tedmirpies,  a  sampling  of  results 
obtained  at  normal  incidence  with  keV  Ax*  ions,  and  others  obtained  at  variable  incident  angles. 
In  additirm,  energy  and  angular  dktributiona  of  ground-state  and  excited-state  sputtered  atoms 
are  compared  in  detail. 

1  Introduction 

Results  of  experimental  studies  of  single  crystals  bombarded  by  keV  particles  offer 
a  unique  perspective  on  the  fundamental  processes  underlying  the  scattering  events 
associated  with  atom  ejection.  Since  Wehner  (19S5)  first  observed  his  anisotropic 
emission  -  his  famous  ‘spots’  -  there  has  been  a  lot  of  spirited  discussion  about 
the  mechanisms  of  «iergy  dissipation,  particularly  with  regard  to  validating  pro¬ 
posed  theories.  While  in  the  early  years  single  crystal  sputtering  was  mainly  of 
interest  for  its  theoretical  significance,  later,  m^e  crystal  effects  became  impor¬ 
tant  in  characterisation  of  surface  structure.  Throughout  this  period,  we  have  seen 
enormous  advances  in  experimental  technique  and  in  the  ability  to  quantitatively 
predict  detailed  information  about  the  trajectories  of  ejected  species.  In  this  chap- 
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in  Isboratflfy  tediBiqaet  which  have  hdpnl  to  reveal  expmmental  evidmce  f<w 
i^tkm  pioceMM  awociated  with  nngie  crystal  targets. 


2  Hbtorical  Perspectives 

The  earliest  systematic  studies  of  the  sputtering  of  single  crystals  date  to  1955  and 
the  pioneering  work  of  Wehner  (1955).  These  studies  showed  that  W  atoms  tgected 
from  nngle-crystal  W  targets  bombarded  by  150  eV  Hg'^'ions  left  the  surface  along 
preferred  crystallographic  directions.  This  anisotropy  was  in  direct  contradiction  to 
many  early  theories  of  sputtering  which  suggested  that  the  process  was  dominated 
by  a  thermal  process  or  by  sublimation  (V<mi  Hippel,  1926  and  Townes,  1944)  which 
would  exhibit  a  yield  with  a  coeine-like  angular  dependence.  Wehner  had  convinced 
his  peers  by  his  use  of  single  crystals  that  momentum  transfer  played  a  key  role  in 
the  phenomenon. 

During  the  1960’s  many  studies  ensued  which  examined  issues  such  as  mass 
dependence,  energy  dependence  and  crystal  face  dep«idence  of  the  sputtering  yidd. 
A  few  key  observations  emerged  from  many  elegant  experiments.  For  the  fee  metals, 
for  example,  it  was  found  that  the  preferred  ejection  generally  occurs  along  the 
<110>  bulk  crystal  direction  (Anderson  U  Wehner,  1960).  Single  crystal  effects 
were  also  observed  for  hep  and  bcc  metals  (Nelson,  1963;  Robinson  k  Southern, 
1968).  Interestingly,  semiconductor  crystals  appeared  much  less  anisotropic  than 
metals  (Southern  et  al.,  1963).  The  reason  for  this  early  erroneous  conclusion 
was  primarily  the  more  extensive  damage  suffered  by  these  materials  from  the 
bombardment  process. 

Single  crystal  sputtering  during  the  early  period  provided  critical  information 
for  the  theoreticians.  As  noted  above,  the  evaporation  models  had  to  be  discarded. 
An  initial  explanation  of  the  anisotropies  came  from  Silsbee  (1957)  who  suggested 
that  momentum  could  be  focused  along  a  row  of  equally  spaced  spheres  if  the 
distance  between  centers  is  less  than  twice  the  atomic  diameter.  This  proposal 
inspired  much  more  thinking  about  the  problem  from  an  atomic  point  of  view.  As 
the  story  goes,  Wehner  attempted  to  look  for  these  Silsbee  chains  by  constructing  a 
model  of  a  crystal  with  balls  attached  by  strings  to  the  ceiling.  He  then  struck  this 
crystal  with  a  primary  particle  (ball)  and  photographed  the  event  with  a  strobe 
lamp.  The  results  were  inconclusive  but  representative  of  bow  clever  Wehner  (1992) 
could  be. 

Computer  simulations  offered  new  hope  to  understanding  single  crystal  sput¬ 
tering.  The  modeling  could  be  performed  on  a  realistic  crystallite  since  the  lattice 
positions  were  known.  Reasonable  interaction  potentials  were  becoming  available 
in  the  mid-1960’8  and  it  was  possible  to  solve  Newton’s  equation  of  motion  for  a 


collection  at(»ns  subject  to  on  initial  particle  bombardment.  Ironically,  some 
of  the  early  modeling  studies  supported  the  conclusions  of  Silsbee  (Gibson  et  al., 
1960).  The  calculations  clearly  revealed  that  focused  collision  sequences  could  am- 
tribute  to  sputtering.  An  early  pioneer  in  computer  modeling,  Don  Harrison,  found 
that  focused  collision  sequences  were  possible,  but  in  fact  were  quite  rare  (Harri¬ 
son  et  al.,  1966,  1968).  More  importantly,  he  suggested  that  the  sputtering  event 
was  dominated  by  near-surface  colliaicms  involving  just  the  first  ccnipk  at  layers 
of  the  solid.  We  will  discuss  the  significance  of  this  important  statement  in  later 
sections.  Lehmann  and  Sigmund  (1966)  later  also  proposed  that  the  surface  struc¬ 
ture  dominated  the  angular  distributkms.  They  noted  that  when  subsurface  atoms 
attain  a  kinetic  energy  on  the  order  of  the  surface  binding  energy,  only  head-on 
coUisiiHU  can  lead  to  emission  of  surface  atoms.  They  also  noted  that  subsurface 
atoms  could  be  qjected  through  potential  minima  in  the  surface  layer,  created  by 
gaps  between  atoms.  These  alternative  explanations  the  angular  anisotropies 
helped  to  explain  a  number  of  problems  with  the  focused  collision  sequence  mod¬ 
els  associated  with  hep  crystals  (Hofer,  1973).  More  importantly,  however,  they 
suggested  that  in  order  to  establish  a  quantitative  understanding  of  smgle-crystal 
sputtering,  it  would  be  essential  for  experimental  data  to  be  obtained  from  crystals 
with  well-characterised  surfaces. 

There  was  much  controversy  about  these  ideas  that  broiled  all  the  way  through 
the  1970’s.  As  knowledge  of  handling  single  crystals  improved,  experimental  meth¬ 
ods  which  required  a  lower  dose  of  incident  ions  were  presented.  For  example, 
Ssymcxak  and  Wittmaack  (1980)  developed  a  sensitive  detection  scheme  for  gold 
using  Rutherford  backscattering  to  analyse  material  collected  on  a  plate.  They 
concluded  that  for  such  high-Z  metals  as  gold  it  might  be  necessary  to  include 
deep-layer  focusing  to  explain  the  large  anisotropies  they  observed.  Hundreds  of 
other  articles  on  single  crystal  sputtering  appeared  during  that  time.  Most  of  these 
have  been  thoroughly  reviewed  by  Hofer  (1991).  He  also  provides  a  comprehensive 
list  of  important  recent  review  articles.  Most  of  this  research  effort  produced  a 
qualitative  understanding  of  single  crystal  sputtering.  For  metals  bombarded  near 
normal  incidence  with  keV  heavy  particles,  the  highest  emission  intensity  is  ob¬ 
served  along  open  crystallographic  directions  in  the  plane  of  the  surface.  Other 
peaks  may  be  observed  in  directions  where  it  is  easy  for  second-layer  particles  to 
escape.  It  was  realized  that  a  high  bombardment  dose  could  produce  topographi¬ 
cal  features  (Hauffe,  1991)  that  could  lead  to  erroneous  results.  Moreover,  surface 
cleanliness  became  an  issue  because  of  the  proposal  by  Harrison,  Lehmann  and 
Sigmund  that  the  sputtering  mechanism  involved  mainly  surface  layers  and  hence 
would  be  influenced  by  adsorbed  impurities.  Energy  distributions  from  single  crys¬ 
tals  were  found  to  be  qualitatively  similar  to  those  found  from  polycrystalline 
surfaces,  although  the  peak  in  the  energy  distribution  may  not  be  the  same  as  for 
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•Bior|>lioiM  qratenw  and  Uim  may  be  atnietuie  ia  tbe  diatributiona  aaaociated  with 
diiect  recoila.  (R«d  at  al.,  1976) 

Mok  quantitative  uaderataBdiog  of  aiogle  ciyatal  aputtering  baa  been  bam- 
peted  by  a  number  of  factors.  First,  experimental  techniques  have  not  been  gen¬ 
erally  available  to  examine  bombarded  single  crystals  at  doses  low  enough  (<10^^ 
iona/cm’)  to  avoid  artifacts.  Secondly,  the  sputtering  community  baa  been  alow 
to  incorporate  auxiliary  surface  characterisation  techniques  such  aa  low  energy 
dectron  dil&action  (LEED)  and  Auger  electron  q>ectroacopy  (AES)  into  their  ap¬ 
paratus.  These  methods  show  that  the  damage  in  single  crystals  is  not  removed 
unless  the  crystals  are  annealed  well  above  their  Debye  temperature  (Jenkins  k 
Chung,  1971).  Moreover,  they  show  that  even  ‘inert’  metals  such  as  copper  and 
gold  cannot  be  kept  clean  unless  the  vacuum  is  reduced  to  the  10~*‘’  tort  regime. 
Unless  single  crystal  surfaces  are  prepared  according  to  the  rigm  of  modem  surface 
physics,  it  is  not  possible  to  be  sure  that  the  resulting  measurements  can  quantita¬ 
tively  be  compared  to  proposed  models  and  that  data  can  even  be  quantitatively 
reproduced  from  lab  to  lab. 

In  this  review,  we  examine  the  presmt  status  of  experimental  approaches  to 
examining  single-crystal  sputtering.  In  general,  the  discussion  will  be  restricted 
to  experimental  results  obtained  on  well-characterised  single  crystal  substrates  ob¬ 
tained  under  low-dose  conditions.  The  emphasis  will  be  to  present  a  limited  data 
base,  primarily  from  our  own  laboratory,  where  quantitative  comparison  to  theory 
is  possible.  It  is  hoped  that  this  review  will  not  only  encourage  more  experiments, 
but  that  the  theory,  now  dominated  by  computationally  intensive  molecular  dy¬ 
namics  calculations,  can  be  reinvigorated  with  new  analytical  approaches. 


3  Experimental  Techniques 

To  fully  characterize  the  sputtering  event  associated  with  single  crystrds,  it  is  desir¬ 
able  to  be  able  to  vary  the  angle  of  incidence  of  the  primary  particle  and  to  detect 
in  a  state-selective  manner  the  sputtered  flux  with  energy,  angle  emd  mass  selection. 
It  is  also  critical  to  obtain  such  data  without  damaging  the  crystal  surface.  During 
the  1970’s  there  were  many  attempts  to  get  at  this  problem.  One  idea  involves  pos- 
tionization  of  the  ejected  species  either  by  electron  bombardment  (Lundquist,  1978) 
or  by  a  low  energy  plasma  discharge  (Oechsner  &  Gerhard,  1972).  The  ion  which  is 
formed  is  then  detectable  within  single  particle  detection  limits  using  a  mass  spec¬ 
trometer.  The  work  of  Thompson  and  his  group  (1963, 1978)  also  stimulated  much 
interest  in  the  time-of-flight  (TOF)  technique.  In  general,  these  methods  require 
an  incident  dose  of  more  than  10^^  particles/cm’,  at  least  3  orders  of  magnitude 
beyond  the  damage  threshold.  Another  approach  involved  the  excitation  of  ejected 
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Fignre  1.  Partial  energy  levd  diagram  for  Rh.  Ground-etate  Rh  can  be  multiphoton  ioniaed 
uaing  two  313.4-nmidiotoiia.  The  fiiat  excited  etate  can  be  probed  by  i<wixing  with  two  338.0-nm 
photoM. 


atoms  by  a  laser,  followed  by  the  measurement  of  the  doppler-shifted  fluorescent 
intensity  (Hammer  et  al.,  1976).  This  scheme  allows  individual  electronic  states 
to  be  detected.  Unfortunately,  the  fluorescence  techniques  have  not  yet  been  sys¬ 
tematically  applied  to  angular  distribution  measurements  and  they  require  beam 
doses  well  in  excess  of  10^^  particles/cm^.  Finally,  an  elegant  experiment  has  been 
devised  to  collect  a  small  amount  of  metal  sputtered  onto  a  detector  using  Ruther¬ 
ford  backscattering  spectroscopy  (Ssymczak  Sc  Wittmaack,  1980).  Here,  the  dose 
could  be  kept  in  the  range  of  10^”  particles/cm^. 

During  the  1980’s  new  laser  techniques  emerged  which  exhibited  the  necessary 
characteristics  needed  for  eflicient  postionisation.  These  include  multiphoton  reso¬ 
nance  ionization  (MPRI)  (Winograd  et  al.,  1982;  Hurst  Sc  Payne,  1988),  multipho¬ 
ton  nonresonance  ionization  (Becker  Sc  Gillen,  1984),  and  single-photon  ionization 
(Schuhle  et  al.,  1988).  The  MPRI  method  seems  best  from  the  point  of  view  of 
maximizing  sensitivity  and  state  selection  as  pointed  out  by  Hurst  et  al.  (1979). 
Other  schemes  may  prove  viable  as  well  in  the  future,  although  none  have  yet  been 
applied  to  the  study  of  the  sputtering  of  single  crystals.  The  MPRI  approach  for 
atoms  can  be  employed  by  inspection  of  the  electronic  energy  level  diagram  shown 
in  Figure  1  for  Rh,  which  will  be  used  as  a  model  system.  There  are  a  variety  of 
electronic  states  that  can  be  directly  probed  by  an  appropriate  choice  of  the  ion¬ 
ization  wavelength.  For  example,  the  ground  state  can  be  ionized  through 
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Figure  3.  Detector  for  perfuming  energy  •i&'l  angle-reeolved  meaetiiemente  of  neutral  epeciee 
deeorbed  from  surfaces,  along  with  the  approp-^’ate  pube  sequences. 


the  *Fr/j  state  with  two  312.4-nin  photons.  The  ^Fy/j,  and  ^Fg/j  states  lie 
within  1  eV  of  the  ground  state  and  can  be  ionized  with  two  photons  of  slightly 
different  energy  (Moore,  1971).  For  this  system  *F9/2*-*Fjf2  in  the  gas  phase  is 
a  forbidden  transition,  so  it  is  possible  to  make  detailed  measurements  on  both 
states  and  to  directly  compare  the  results. 

The  experimental  geometry  for  such  a  scheme  is  shown  in  Figure  2  (Winograd 
et  al,  1986).  A  pulsed  ion  beam  of  200  ns  serves  to  generate  a  burst  of  desorbed 
particles  localized  in  space  moving  away  from  the  target  surface.  At  an  ^propriate 
time,  a  laser  pulse  tuned  to  an  atomic  resonance  of  the  emitted  particle  is  directed 
a  few  millimeters  above  and  parallel  to  the  target.  As  discussed  by  others  (Hurst 


k  Payne,  1988),  with  nifficient  laser  power  and  excited-state  lifetimes,  all  neutral 
spedes  with  excitations  in  resonance  with  the  photon  field  and  with  appropriate 
ionisatkm  potentials  are  converted  quantitatively  into  ions.  Due  to  the  pulsed 
natate  the  expoiment  and  the  selective  nature  of  the  ionization  process  itself, 
the  resulting  ions  may  be  efilciently  counted  by  accelerating  them  with  an  external 
potential  and  measuring  their  TOF. 

In  this  configuration  the  laser  beam  is  focused  in  the  shiq>e  of  a  ribbon  0.1 
mm  thick  and  6  mm  wide  and  directed  1  cm  above  and  parallel  to  the  target. 
Focunng  allows  kuiisation  of  those  neutral  species  which  are  in  a  specific  region 
of  apace.  Hence,  by  varying  the  time  difference  between  a  200-ns  ion  pulse  and 
a  fi-ns  laser  pulse,  the  velocity  distribution  of  the  neutral  species  is  determined. 
Fhrthermore,  the  resulting  ions  may  be  imaged  onto  a  microchannel  plate  which 
allows  the  pontion  of  the  neutral  species  in  the  laser  beam  to  be  determined  with 
respect  to  the  pontion  of  the  incident  ion  beam  on  the  target.  With  knowledge  of 
the  nature  of  the  extraction  fields  and  instrumental  factors  associated  with  TOF 
measurements  and  position-sensitive  detectors,  it  is  feasible  to  calculate  the  take-off 
angles  of  the  neutral  species  (Kobrin  et  al.,  1986).  This  iq>paratus  was  also  suitably 
equipped  with  LEED  for  surface  crystallography  measurements  and  AES  for  surface 
cleanliness  monitoring.  This  energy-  and  angle-resolved  neutral  (EARN)  detector 
is  a  powerful  device  for  measuring  neutral  trajectories  with  enormous  sensitivity. 
Typically,  a  full  EARN  measurement  can  be  recorded  with  a  dose  of  10^^  Ar*** 
ioBs/cm^. 

4  Experimental  Studies  of  Single  Crystal  Bom¬ 
bardment 

4.1  Yield  Measurements 

Although  extensive  effort  has  been  devoted  to  measuring  the  total  yield  of  particles 
emitted  from  polycrystalline  targets  bombarded  at  normal  incidence,  not  many  of 
these  measurements  have  been  reported  on  single  crystals  since  the  1960’s  (Mag- 
nuson  k  Carlston,  1963).  These  workers  used  a  weight  loss  technique  to  determine 
that  the  yield  of  Cu{100}  is  2.0  times  the  yield  of  Cu{110}  and  0.7  times  the 
yield  of  Cu{lll}.  This  trend,  generally  supported  by  computer  simulations  (Har- 
risoD  et  al.,  1979),  follows  the  trend  of  surface  atom  density  in  these  materials.  In 
general,  these  measurements  are  quite  variable  due  to  the  ill-defined  experimental 
conditions  and  are  of  only  qualitative  value. 

An  important  recent  measurement  using  nonresonant  laser  postionization  pow¬ 
erfully  illustrates  the  necessity  of  using  clean  well-defined  surfaces  under  low  dose 
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Figure  3.  The  Ru  SNMS  eignul  ratio,  —  0)>  from  the  Ru(OOOl)  mirface  plotted  ae 

a  ftmction  of  primary-ioa  doee,  J.  Reeult*  from  four  leparate  experiments  are  shown  in  (a)  while 
results  from  a  single  experiment  are  shown  in  (b).  The  solid  lines  are  the  result  of  least-squares 
fits  of  the  data  by  a  proposed  model  (FVom  Burnett  et  al.,  1989). 

conditions.  As  shown  in  Figure  3  the  yield  from  Ru(OOOl)  decreases  by  a  factor 
of  two  as  the  dose  is  increased  above  ~  10^^  ions/cm^  for  bombardment  with  3 
keV  Ar***  at  normal  incidence  (Burnett  et  al.,  1989).  Since  most  other  techniques 
require  at  least  1x10^^  particles/cm^,  it  is  clear  that  those  studies  are  tainted  by 
significant  surface  damage  effects.  For  the  Ru  case,  it  is  tentatively  proposed  that 
the  decrease  in  yield  arises  from  the  primary  ion  going  into  defect  sites  created  by 
previous  bombardment.  As  seen  in  the  Figure,  this  model  fits  well  with  the  ex¬ 
perimental  data,  although  there  are  a  number  of  adjustable  parameters  in  it.  This 
idea  needs  to  be  further  tested  since  it  does  not  seem  likely,  based  on  the  results 
of  recent  molecular  dynamics  computer  simulations  (Harrison,  1988).  It  would  be 
valuable  to  see  similar  studies  on  other  materials  and  at  various  temperatures. 

4.2  Depth  of  Origin  of  Emitted  Particles  from  Single  Crys¬ 
tals 

In  order  to  have  a  proper  mechanistic  view  of  the  sputtering  event,  it  is  of  course 
essential  to  be  certain  about  from  where  the  emitted  particles  originate.  Contro- 
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Figure  4.  Schemutk  repreeentetion  of  our  edieiiie  far  repreecnfing  ciystel  direction*  far  on 
fec{lll}  and  {001}  face.  A  polar  aagk  d  s  0*  k  normal  to  tlie  eurface.  Open  circle*  deaig- 
nate  first-  l^rer  atom*  and  ahaded  cirde*  aecond-lajwr  atom*.  The  letter*  A-C  dekgnate  po**ible 
adaorption  nte*  for  oxygen  atom*. 


veny  about  this  subject  dates  to  the  original  model  by  Sigmund  (1969)  where  it 
was  implied  that  particles  could  be  emitted  from  below  the  surface  layers.  Later, 
computer  simulations  (Harrison  et  al.,  1978)  showed  that  for  copper  single  crystals 
>90%  of  the  particles  are  emitted  from  the  top  layer.  A  corrected  sputter  depth 
has  recently  been  given  (Vtcanek  et  al.,  1989)  which  is  conristent  with  simulations. 

Experimental  tests  of  this  idea  on  monoelemental  targets  are  not  available. 
Ideally,  it  would  be  interesting  to  create  a  single  crystal  target  with  an  isotope 
whose  surface  layer  is  different  from  the  bulk  and  study  it  using  laser  postionization 
detection.  Several  workers  have  attempted  to  come  close  to  accomplishing  this  goal. 
Prigge  and  Bauer  (1980)  measured  the  W**  yield  from  W{  110}  covered  with  varying 
amounts  of  Cu,  Ag  and  Pd.  Although  the  results  are  subject  to  many  artifacts, 
their  studies  show  that  the  W***  signal  disappears  between  one  and  two  monolayers 
of  deposited  overlayers. 

The  most  definitive  experiment  on  well-defined  single  crystals  was  carried  out  on 
Ru{0001}  covered  with  known  amounts  of  Cu  (Burnett  et  al.,  1988)  and  bombarded 
at  normal  incidence  by  3.6  keV  Ar'*'  ions.  The  experiment  is  particularly  nice  since 
they  employed  low  dose,  LEED/Auger  and  used  MPRI  postionization  for  detection. 
Moreover,  Cu  is  known  to  form  two-dimensional  rafts  on  Ru,  rather  than  multilayer 
structures.  The  results  clearly  show  that  virtually  all  of  the  emitted  particles  arise 
from  the  first  layer. 

Is  this  the  final  answer?  No,  there  are  a  number  of  systems  where  the  depth  of 
origin  is  probably  more  than  one  layer.  Computer  modeling  of  Si{001 } ,  for  example, 
(Smith  et  al.,  1989)  clearly  shows  that  particles  can  be  emitted  from  the  fourth  or 
fifth  layer.  There  is  no  direct  experimental  evidence  for  this  prediction  yet,  but 
recent  studies  on  GaAs{001}  by  Burnham  et  al.  (1993)  indicate  that  fourth  layer 
Ga  atoms  can  contribute  significantly  to  observed  angular  distributions.  Recent 
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Figure  5.  EiARN  intcnaity  maiM  for  dean  IUi{lll}  and  p(2x2)  O  IU>{111}.  The  ploU  arc 
nonBalned  to  the  highest  intensity  peak  in  both  cases.  The  podtive  values  of  #  ate  recorded 
along  d  —  +30*  and  the  negative  values  of  t  are  recorded  along  4  ^  —30*  (FVom  Winograd  ct  al., 
19M).  Note  that  in  the  older  spot  notatioa,  the  [110]  spot  would  occur  at  9  s  35.3*  an^  s  -30* 
and  the  [100]  spot  would  occur  at  d  s  54.7*  and  4  —  +30*. 


studies  using  the  scanning  tunneling  microscope  (STM)  also  suggest  that  deeper 
layers  are  exposed  by  single  ion  impacts  on  Si  single  crystals  (Feil  et  al.,  1992).  If 
these  results  are  correct,  they  may  be  important  in  understanding  depth  profiling 
in  these  commercially  important  materials. 


4.3  Energy  and  Angular  Distributions  from  Metallic  Single 
Crystals 

As  was  true  during  the  time  of  Wehna’s  early  experiments,  accurate  EARN  dis¬ 
tributions  from  well-defined  surfaces  are  important  for  testing  various  theoretical 
models.  Rhodium  represents  an  important  prototypical  system  for  model  studies. 
As  seen  in  Figure  1,  the  electronic  structure  of  Rh  is  well-suited  for  state-selected 
MPRI  studies.  Moreover,  the  surface  chemistry  of  Rh  has  been  extensively  investi¬ 
gated  by  many  techniques,  and  the  surfaces  are  known  to  exhibit  minimal  surface 
relaxation  or  reconstruction  that  might  complicate  the  analysis  (Shepherd  et  al., 
1978). 

Let  us  first  consider  the  (111)  face  bombarded  at  normal  incidence  by  5  keV 
Ar***  ions  (Winograd  et  al.,  1986).  The  crystallographic  definitions  for  this  face 
are  defined  in  Figure  4,  and  the  EARN  data  are  shown  in  Figure  5.  Note  that  for 
this  crystal  face,  the  surface  atoms  exhibit  six-fold  symmetry.  The  three  second- 
layer  atoms,  however,  create  bulk  three-fold  symmetry  with  characteristic  direc¬ 
tions  along  <111>  (^  =  0*  in  our  notation),  <211>  4  =  -30*,  and  <112>, 
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Figure  6.  Polar  angle  distributions  for  various  asimuthai  angles  for  fixed  secondary  kinetic  energy 
of  the  Rh  atoms.  In  each  frame  the  data  are  normalised  to  the  ^  =  -30*  peak  intensity.  The 
calculated  data  using  the  EAM  potential  are  reported  with  a  full  width  at  half-maximum  (fwhm) 
resolutitmof  IS*  in  the  polar  angle.  A  constant  solid  ani^e  is  used  in  the  histogramimng procedure. 
The  experimental  resolution  if  also  approximatriy  IS*.  The  surface  normal  corresponds  to  S  s 
0*. 


^  =  -f-30*.  Experimental  EARN  results  for  ^  ±  30**  are  shown  in  Figure  5.  Note 
that  the  results  are  qualitatively  similar  to  many  of  the  early  data  taken  under 
much  more  trying  circumstances.  In  Wehner’s  nomenclature,  there  would  be  a  set 
of  3  weaker  spots  =  30”,  150*  and  270*),  3  stronger  spots  =  -30*,  90*  and 
210*),  as  well  as  a  very  weak  central  spot.  The  intensity  along  the  closest-  packed 
line  of  atoms  in  the  surface  plane  (4  =  0*,  60*,  ■  •  •)  is  not  shown  but  exhibits 
a  minimum  in  intensity.  The  energy  distributions  are  clearly  dependent  on  the 
take-off  angle,  suggesting  that  their  final  energies  are  influenced  by  surface-layer 
collisions. 

This  directional  nomenclature  is  different  from  that  previously  associated  with 
single  crystal  sputtering.  The  reason  for  this  change  is  that  the  older  notation  was 
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bawd  on  the  aasampti<M)  that  the  uiaotropy  aroee  primarily  from  focuaed  coUiricm 
sequences  deep  within  the  crystal.  Hierdbre,  it  was  natural  to  define  the  ‘spots’ 
using  bulk  crystallographic  directions.  Now,  we  know  that  surface  processes  are 
the  main  contributors  to  the  angular  anisotropies.  It  is  essential,  therefore,  that 
we  adopt  a  notation  that  reflects  this  process  and  is  consistent  with  the  notation 
employed  by  the  surface  physics  community. 

With  present-day  computers,  it  is  possible  to  utilise  molecular  dynamics  com¬ 
puter  simulations  to  try  to  reproduce  the  experimental  results  shown  in  Figure 
5.  These  calculations  are  important  to  carry  out  for  two  reasons.  First,  they  tell 
us  whether  a  molecular  dynamics  model  is  sufflcient  to  describe  the  formation  of 
the  collision  cascade  and  the  subsequent  trajectories  of  ejecting  atoms.  Second, 
they  allow  us  to  examine  the  microscopic  mechanisms  that  give  rise  to  the  angular 
anisotropies.  The  comparison  between  theory  and  experiment  provides  a  modern- 
day  reality  check. 

The  details  of  the  molecular  dynamics  scheme  in  general  are  described  elsewhere 
in  this  book  (Robinson,  1993;  Nieminen,  1993).  The  details  of  the  procedure  for 
calculating  EARN  distributions  for  Rh  have  been  described  in  detail  and  will  not 
be  discussed  here  (Garrison  et  al.,  1988;  Foiles  et  al.,  1986).  As  seen  in  Figure  6, 
however,  it  is  clear  that  for  Rh{lll},  at  least,  the  comparison  between  calculation 
and  experiment  is  about  as  good  as  can  be  expected  over  a  wide  rang''  of  phase 
space. 

The  calculations  tell  us  that  the  collision  sequences  are  strongly  dominated  by 
the  alignment  of  atomic  motions  inside  the  solid.  As  these  motions  cause  ejection  of 
first-layer  atoms,  further  focusing  is  caused  by  channeling  or  blocking  by  other  first- 
layer  atoms.  Note  that  this  type  of  channeling  occurs  at  lower  energies  than  that 
often  associated  with  higher  speed  phenomena.  For  example,  the  highest  intensity 
is  observed  along  the  open  crystallographic  directions  —  ±  30*  in  our  case),  and 
the  minimum  intensity  is  observed  along  the  close-packed  crystallographic  direction 
(0=  0*).  If  only  the  top  surface  layer  was  important,  the  peaks  at  9  =  -37*,  ^  =: 
-30*  and  6  =  +42* ,  ^  =  +30*  should  be  equal  in  intensity  and  not  unequal  as  shown 
in  Figure  5.  The  additional  intensity  at  0  =  -37*,  ^  =  -30*  arises  mainly  from  the 
ejection  of  atom  A  by  atom  B  (Figure  4)  with  first-layer  focusing  by  other  surface 
atoms.  The  peak  at  0  =  +42*  and  ^  +30*  is  lower  in  intensity  by  a  factor  of 

0.5  at  low  kinetic  energy  since  no  such  mechanism  is  available  along  this  azimuth. 
The  peak  at  9  =  0*  has  a  large  component  from  ejection  of  the  second-layer  atom 
B  which  is  focused  upward  by  three  surface  atoms. 

Many  earlier  experiments  have  been  reported  on  {111}  surfaces.  Szymczak  Si 
Wittmaack  (1980)  found  for  4  keV  Ne***  on  gold,  for  example,  that  the  peak  along 
^  =  -30*  (the  [Oil]  spot  in  their  notation)  was  at  least  4  times  the  intensity  of  the 
peak  along  4>  =  +30*  (the  [100]  spot).  No  intensity  above  a  cosine  background  was 
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obwrved  at  normal  emission.  Hie  results  are  important  since  they  show  that  the 
overall  intensity  directed  into  the  qiots  is  quite  hi|^,  up  to  50%  in  favorable  cases. 
They  are  quantitatively  different  from  the  results  found  for  Rh{lll},  however,  in 
that  there  is  less  of  a  difference  found  between  4  =  -30*  and  ^  =  +30*  and  there 
is  a  small  amount  of  intensity  found  at  normal  emission  for  Rh{lll}.  It  is  unlikely 
that  the  distributions  reported  by  Ssymcsak  and  Wittmaack  could  be  reproduced 
by  the  computer  modeling  as  these  calculations  would  probably  yield  results  very 
siinilar  to  those  found  on  Rh{lll}.  Without  further  experiments  and  calculations, 
however,  further  hypothesising  about  these  differences  would  be  mere  speculation. 

Another  study  of  Ni{lll}  bombarded  by  1.1  keV  Ne*^  ions  at  normal  incidence 
has  recently  been  reported  by  Wucher  et  al.  (1992).  Their  angular  distributicms, 
which  are  energy  selected,  are  very  similar  to  those  found  on  Rh{lll}.  For  exam¬ 
ple,  they  find  that  only  the  relatively  high  energy  particles  are  emitted  at  normal 
incidence  as  shown  in  Figure  6.  The  major  difference  arises  from  the  observation 
that  the  A  =  +30*  peak  is  sevml  times  more  intense  than  the  ^  =  -30*  peak  after 
the  background  is  removed.  Th^  computer  modeling  suggests  these  two  features 
should  be  much  closer  in  intensity  as  suggested  by  the  data  for  Ilh{lll}  shown  in 
Figure  6.  Since  their  experiments  were  not  performed  in  the  low-dose  regime,  it  is 
difficult  to  speculate  about  the  cause  of  this  discrepancy. 

How  do  these  measurements  and  interpretations  jive  with  earlier  models  of  the 
production  angular  anisotropies?  This  is  a  very  tricky  question.  Mechanis¬ 
tic  conclusions  drawn  from  molecular  dynamics  calculations  involve  an  analysis  of 
thousands  of  trajectories  in  an  attempt  to  extract  the  essential  physics.  As  Har¬ 
rison  has  often  illustrated  (1988),  the  atom  ejection  process  is  associated  with  a 
cataclysmic  event.  Although  explanations  of  the  angular  anisotropies  for  Ilh{lll} 
arise  from  such  an  analysis,  there  are  many  other  mechanisms  which  also  contribute 
significantly  to  sputtering.  It  is  often  not  possible  to  distill  all  this  particle  motion 
into  simple  pictures.  The  critical  issue  here  is  that  the  molecular  dynamics  cal¬ 
culations  can,  for  the  first  time,  provide  quantitative  agreement  with  experiment, 
at  least  for  Rh{lll}.  The  important  mechanisms  are  embedded  in  the  computer 
printouts  and  are  available  for  the  taking  even  though  with  considerable  difficulty. 
This  degree  of  certainty  has  never  been  available  before. 

Given  this  caveat,  there  appears  to  be  a  great  deal  of  similarity  between  the 
notions  originally  presented  by  Harrison  et  al.  (1966),  by  Lehmann  and  Sigmund 
(1966)  and  by  the  most  recent  studies  -  as  long  as  you  do  not  look  too  closely. 
The  early  work  focused  attention  on  mechanisms  involving  transparency  effects, 
lens  effects  and  specific  force  fields.  Both  Lehmann  and  Sigmund  mecbruiisms  are 
quite  reasonable  in  light  of  the  molecular  dynamics  studies,  although  they  begin  to 
break  down  when  you  ask  questions  about  the  energy  dependence  of  the  angular 
distributions.  Their  model  would  also  predict  that  there  should  be  no  intensity 


226 


lifM4S 


Figure  7.  CroM-MctkMiel  view  of  the  Rh{331}  nufece. 


in  the  ^  =  +30*  direction.  The  key  question,  then,  involves  the  generality  of  the 
mechanisms  discussed  above  when  applied  to  a  variety  of  other  systems. 

If  channeling  and  blocking  in  the  first  layer  ue  important  in  controlling  the 
angular  distributions,  there  should  be  similar  trends  found  on  other  crystal  faces. 
Rh{331}  is  an  interesting  surface  since  it  is  a  stepped  structure  with  a  (111)  terrace 
that  is  three  atoms  wide.  It  is  very  similar  to  the  {111}  face  except  that  it  is  tilted 
by  21.9*  as  shown  schematically  in  Figure  7.  The  experimental  results  show  that 
the  ejection  is  strongly  peaked  along  the  4  =  +90*  asimuth,  and  that  desorption 
along  other  crystallographic  directions  is  consid^ably  reduced  (Reimann  et  al., 
1989).  Molecular  dynamics  calculations  show  that  the  9  =  15*  peak  arises  from 
the  same  surface  channeling  mechanism  operative  for  {111}  corresponding  to  the 
0  =  37*  peak  along  ^  =  0*.  The  reduced  intensity  along  the  other  directions  arises 
either  from  the  presence  of  open  channeb  in  the  crystal  (near  9  s  0*)  or  from 
blocking  due  to  the  atomic  step  along  4  =  -90*. 

Similar  anisotropies  are  observed  for  Rh  atoms  desorbed  from  ion>bombarded 
Rh{001}  (Maboudian  et  al.,  1990;  El-Maazawi  et  al.,  1991).  The  experimental 
EARN  distribution  for  this  surface  is  shown  in  Figure  8.  The  highest  intensity  is 
found  along  the  open  crystallographic  direction,  ^  =  0*,  while  a  minimum  intensity 
is  observed  along  4  =  45* .  A  comparison  between  the  polar  angle  distributions  for 
Rh{lll}  and  Rh{001}  reveals  additional  details  about  the  surface  channeling  and 
blocking  mechanisms.  It  turns  out  that  the  peak  in  the  polar  angle  distribution 
in  the  10-20  eV  range  is  larger  for  {001}  (dmmx  =  40*,  4  =  0*)  than  for  {111} 
(^max  =  37*,  4  =  -30*).  This  difference  arises  since  the  open  channel  is  larger  on 
{001},  allowing  more  grazing  take-off  angles.  The  polar  angle  intensity  maxima 
along  the  close-packed  directions,  however,  are  nearly  identical  for  both  crystal 
surfaces  (^mwc  =  34*  for  {001}  and  9max  =  32*  for  {111}).  This  result  is  important 
since  the  distance  to  the  nearest-neighbor  atom  along  the  close-packed  direction 
is  identical  for  both  surfaces.  Thus,  the  polar  angle  distributions  are  determined 
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Fifore  8.  Enargy-  and  an^»-rc«alved  dktributioM  of  Rh  atooM  in  Um  *Ft/2  axdted  rtata  and 
*F2/2  ground  atata,  qoctod  from  5  keV  Ar'*'  ion  bombardod  fUi{100}.  The  data  comqMnd  to 
ajedion  along  ^  s  0*  (<100>)  and  m  45*  (<110>)  ciyataIlo(p^>hic  dnectiona,  aa  ddBned  in 
the  inaet.  Due  to  the  qrnunetiy  of  the  auiiace  and  the  angular  readution  (e.g.,  s  dk  8*  at  8 
s  45*),  the  reauha  repaeaent  the  data  over  ~  50%  of  aU  qMce.  Both  plota  are  nonnaliied  to  the 
maximum  intenaitf  peaka.  See  Figure  4  for  definitiona. 


by  the  extremely  local  interactions  experienced  by  the  desorbing  atom,  and  can  be 
predicted  from  the  presence  or  absence  of  surface  channeling  directions. 

Since  crystal  structure  via  surface  channeling  and  blocking  strongly  influences 
the  EARN  distributions,  it  follows  that  adsorbate  atoms  or  molecules  on  single¬ 
crystal  substrates  should  systematically  alter  the  trajectories  of  the  desorbing  un¬ 
derlayer  species.  The  ctmcept  has  been  tested  in  detail  for  the  p(2x2)  O  ordered 
overlayer  on  Rh{lll}  with  the  EARN  experiment  (Winograd  et  al.,  1986;  Reimann 
et  al.,  1989).  For  an  atomic  adsorbate  of  this  sort,  there  are  a  number  of  possible 
high-symmetry  binding  sites  including  two  different  3-fold  hollow  sites  (often  re¬ 
ferred  to  as  the  B  or  hep  site  and  the  C  or  fee  site)  and  an  on-top  site  (or  A  site) 
as  diown  in  Figure  4.  The  simplest  idea  is  that  adsorption  of  oxygen  atoms  at  a  B 
site  should  preferentially  block  Rh  atoms  desorbing  along  the  ^  =  -<-30*  asimuth, 
adsorption  at  a  C  site  should  alter  Rh  atoms  leaving  along  ^  =  -30°  while  A  site 
adsorption  may  have  very  little  effect,  except  perhaps  on  the  particles  emitted  at  9 
s  0*.  Polar  an^e  measurements  as  seen  in  the  lower  part  of  Figure  4  clearly  show 
that  the  ^  =  -30*  direction  is  preferentially  reduced  relative  to  either  9  =  0*  or  ^ 
s  -<-30*  stroni^  suggesting  C  site  adsorption. 

These  results  are  quite  interesting  since  they  clearly  show  how  sensitive  the 
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imoifition  anfolar  dtstiibatioBt  we  to  mieU  ovorlayer  coverages.  In  this  caae, 
in  adiUtion  to  the  obaecved  pertarbatkMM  in  iht  EARN  dwtcibutiona,  the  total 
ground  state  neutral  Rh  atom  yield  was  observed  to  decrease  by  a  factor  of  about 
2  (Reimann  et  al.,  1989).  These  observations  further  illustrate  that  tri^M^y 
measurements  which  are  to  be  cconpared  with  theory  must  be  performed  on  well- 
cleaned  and  characterised  surfaces  and  under  low-dose  conditions. 

4.4  Energy  and  Angular  Distributions  from  Semiconductor 
Single  Crystals 

Angular  anisotropies  have  also  been  observed  for  atoms  qected  from  ion-bombarded 
single  crystal  semiconductor  surfaces.  The  earliest  observations  (Southern  et  al., 
1963)  indicated  that  the  anisotropies  were  smaller  than  for  metals.  Later,  however, 
it  was  discovered  that  the  angular  disttibutions  were  strongly  temperature  depen¬ 
dent,  suggesting  that  crystal  damage  and  thermal  annealing  were  complicating 
the  analysis  (Anderson  et  al.,  1963).  Because  of  the  technological  importance  of  Si, 
GaAs  and  related  materials,  many  studies  of  the  topological  changes  that  occur  due 
to  extensive  bombardmoit  have  been  reported  (Carter  et  al.,  1984).  Only  recently 
have  well-defined  surfaces  prepared  under  UHV  conditions  been  studied  in  detail. 
What  emerges  is  a  fascinating  comparison  between  the  dynamics  of  sputtering  in 
metals  with  that  found  in  the  more  open  covalent  lattice  of  semiconductors. 

GaA8{001}  represents  an  extraordinary  model  semiconductor  surface.  Depend¬ 
ing  upon  how  it  is  grown,  the  surface  can  be  terminated  with  various  coverages  of 
As  or  Ga  (Arthur,  1974;  Neave  k  Joyce,  1978).  It  crystallises  with  a  zinc  blende 
lattice,  the  same  as  Si,  which  means  that  the  001  plane  exhibits  alternating  layers 
of  As  and  Ga  atoms.  Moreover,  the  surfaces  exhibit  a  well-defined  rearrangement 
that  has  been  characterized  by  scanning  tunneling  microscopy  (Pashley  et  al.,  1988) 
and  by  various  diffraction  techniques.  The  surfaces  are  normally  prepared  using 
molecular  beam  epitaxial  growt^methodology  (Arthur,  1974).  ^  KB  B) 

Molecular  dynamics  computer  simulations  have  only  recently  been  attempted 
for  semiconducting  materials  (Stansfield  et  al.,  1989;  Smith  et  al.,  1989).  The 
calculations  are  more  difficult  in  the  sense  that  there  is  directional  bonding  which 
requires  the  use  of  many-body  potential  functions.  On  the  other  hand,  the  lattices 
are  much  more  transparent  to  the  primary  particle  so  that  nuclear  positions  play  a 
more  dominant  part  in  controlling  the  dynamics  than  the  chemical  bonding  forces. 
The  fact  that  the  knowledge  of  the  dynamics  is  not  as  extensive  for  semiconductors 
as  it  is  for  metals  provides  an  added  incentive  for  more  effort  in  this  area. 

It  is  easy  to  examine  the  influence  of  the  structure  of  the  first  atomic  layer 
on  the  angular  distribution  of  sputtered  atoms  from  inspection  of  Figure  9.  This 
figure  shows  on  the  left  side  a  hypothetical  unreconstructed  GaAs{001}  arsenic- 
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Figure  9.  T<q>  paael:  RepwieentHioni  of  three  hypothetical  G*As{001}  euifacea.  The  yellow  ball* 
repreaent  eurtaoe  Aa  atone.  The  purple  baUs  repreaent  aecond-layer  Qa  atoma  and  the  green  balla 
repreaent  third  layer  Aa  atonu.  Bottom  panel:  calculated  angular  diatiibuticma  of  10>30h!V  Ga'*’ 
iona  deaorbed  by  keV  Ar'*’>ion  bombardment  of  the  ccureaponding  GaAa{001}  aurface  given  in 
the  top  panel.  Ibe  <011>  aaimuthal  cryatal  direction  ia  parallel  to  arrow  (c)  in  the  figure.  The 
polar  angle  ^  Ga'^-ion  emiaaion  ia  proportional  to  the  diatance  a  q>ot  from  the  center  of  the 
circle.  Ion  oiectioB  mechaniama  (a),  (b)  and  (c)  are  diacuaaed  in  the  text.  (Fh»i  BInmenthal  et 
al.,  1991). 


terminated  surface  consistuig  of  a  square  array  of  As  atoms  bonded  to  two  Ga  atoms 
in  the  li^er  below.  Each  As  atom  possesses  two  partially  filled  ‘dangling  bonds’ 
pointing  upwards  and  oriented  parallel  to  the  <001  >  directions.  The  As  atoms 
may  pair  up  to  form  dimers  along  this  direction,  doubling  the  lattice  periodicity, 
also  shown  in  Figure  9.  In  the  laboratory,  the  (2x  1)  reconstruction  is  not  observed 
f<»  GaAs  but  is  for  Si  and  C.  Instead  a  series  of  structures  is  found,  the  most 
important  of  which  is  the  As-rich  (2x4)  geometry  shown  in  the  right  most  panel. 
This  structure  is  found  during  MBE  growth  of  single-crystal  GaAs  films  (Pashley 
et  al.,  1988;  Biegelsen  et  al.,  1990). 


imf  4s 


Figure  10.  Top;  The  direct  gectioamochoniMnwIiicfaiefoepeneible  for  the  two  Auniniit  toot ure« 
ia  the  GoAe{001}-(2x4)  eecoodory  ion  diatributioo.  The  meehaainn  it  dtetacteiiaed  by  a  direct 
moeneatum  tranefer  between  a  third  layer  Ae  atom  and  a  aecoad  layer  Ga  atom  aloag  their  bond 
axie  in  the  <110>  direction.  The  Ga  atom  may  qect  over  the  ehaand  between  Aaa  dimen  or 
up  through  the  Aei  dimen  aa  pictured.  Bottmn:  The  deeorption  of  a  eecond  l^rer  Ga  atom  in 
the  <110>  direction  acroee  the  mieeing  dimer  row.  The  pictured  coUieion  sequence  ehowe  a  lower 
layer  Ae  atom  channeling  up  throu^  the  crystal  and  tranafetring  momentum  to  the  Ga  sttom. 
The  Ga  atom  may  desorb  at  a  high  polar  angle  since  it  ejects  across  the  missing  dimer  row.  The 
momentum  transfer  in  this  direction  is  not  as  eflSdent  as  in  the  <110>  direction  which  causes 
the  desorbed  ion  intensity  to  be  lower. 


Molecular  dynamics  calculations  are  not  yet  amiable  for  GaAs.  It  is  possible, 
however,  to  approximate  the  GaAs  response  to  bombardment  with  a  Si  lattice  since 
they  both  have  related  bonding  properties  and  identicd  crystal  structures.  Results 
of  molecular  dynamics  calculations  for  Si{001},  Si{001}-(2xl)  and  Si{000}-(2x4) 
are  shown  in  the  bottom  part  of  Figure  9.  The  computational  details  have  been 
discussed  elsewhere  (Smith  et  al.,  1989).  The  distributions  are  displayed  in  a  format 
which  allows  a  real-space  comparison  to  the  model  crystals  of  GaAs  shown  above 
them.  The  asimuthal  angles  are  defined  relative  to  the  center  of  each  circle.  The 
polar  angle  is  related  to  the  distance  of  a  spot  from  the  center  of  the  circle.  These 
distributions  are  for  the  second-layer  Si  atoms  which  are  found  to  eject.  This  set 
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Flgurt  11.  TIm  rtUtive  intaaaity  of  30-cV  Q*'^  ion  dotorbod  by  S-keV  nonnal  incident  Ar'^-ion 
bombardment  of  Ute  QaAe{001}>(3)(4)  euilace.  The  polar  anfle  ie  the  angle  of  detection  from 
the  mirface  normal.  For  the  in-plmte  aaimuthal  angle,  4  —  0*  correepondr  to  the  <011>  ctyetal 
dfawetioa. 


ig  crystnllogr^hicadly  equivalent  to  the  Ga  atoms  shown  by  the  shaded  atoms  of 
Figure  1. 

In  a  qualitative  sense,  the  results  of  these  calculations  are  quite  fascinating. 
For  the  unreconstructed  GaAa{001}  surface,  the  Ga  distribution  exhibits  only  two 
peaks,  at  ^  =  90*  and  270*,  with  no  structure  found  at  higher  polar  angles  and 
no  intensity  at  ^  s  0*  and  180*  or  9  <45*.  A  detailed  analysis  of  the  atomic 
trajectories  leading  to  these  distributions  suggests  that  the  anisotropy  arises  from 
a  direct  mechanism  whereby  a  third-layer  As  atom  collides  with  a  second-layer 
Ga  atom  and  thereby  causes  the  Ga  atom  to  eject  along  their  mutual  bond  axis 
(mechanism  a).  This  mechanism  can  be  ascribed  to  the  directional  bonding  and 
open  structure  of  covalent  crystals  and  is  rarely  observed  on  metal  surfaces.  A 
graphical  example  of  this  mechanism  is  shown  in  Figure  10  at  the  top.  The  effect 
of  the  (2x1)  and  (2x4)  reconstructions  is  also  clearly  visible  through  a  broadening 
of  the  major  spots  (mechanism  b)  and  by  the  appearance  of  intensity  at  ^  =  0° 
and  ^  =  180°,  respectively.  The  mechanisms  associated  with  these  processes  are 
somewhat  varied,  but  have  all  been  identified.  A  graphical  example  of  mechanism 
c  is  shown  in  the  bottom  of  Figure  10. 

The  experimental  results  for  Ga'*'  ions  sputtered  from  GaAs{001}-(2x4)  bom¬ 
barded  by  1  keV  Ar'*'  ions  are  shown  in  Figure  11.  Although  the  laser  postionisation 
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Figure  12.  Left:  Tbe  tr^jcctoiiei  bom  tbe  inteructioB  of  e  3  keV  Ar'^  ion  beam  with  a  Ni  atom. 
The  inckkiit  iona  are  deflected  and  a  legioii  ol  aero  ion  flux  ie  created  behind  the  Ni  atom.  The 
edge  of  thia  region  definea  the  ahadow-cone.  All  diatancea  are  in  Angatidma.  Right:  A  profile  of 
the  ahadow>coiie  5  A  bdiind  the  Ni  atom  which  illuatratea  the  focuaing  eflect. 


experiments  have  yet  to  be  performed  on  this  surface,  related  studies  on  GaAs{  110} 
suggest  that  the  ion  and  neutral  distributions  are  quite  similar  (Blumenthal  et  al., 
1990).  Note  that  the  general  features  predicted  by  the  calculations  are  clearly 
visible  and  that  the  influence  of  the  reconstruction  can  be  seen  at  higher  polar  an¬ 
gles  as  expected  (Blumenthal  et  al.,  1991).  These  experiments  are  continuing  with 
a  special  emphasis  on  how  the  features  should  change  during  deposition  of  thin 
metallic  overlayers  such  as  Al.  It  would  also  be  interesting  to  be  able  to  measure 
the  As  emission,  although  at  the  present  time  it  is  experimentally  intractable. 


5  Angle  of  Incidence  Effects  -  Shadow- Cone  En¬ 
hanced  Desorption 

Much  of  what  we  have  learned  about  angular  distributions  so  far  is  derived  from 
experiments  carried  out  at  normal  incidence,  partly  because  the  dynamics  cal¬ 
culations  are  easier  under  these  conditions,  and  partly  because  of  experimental 
convenience.  It  is  well-known  from  the  early  literature,  however,  that  tbe  angle  of 
incidence  of  the  primary  particle  significantly  influences  the  sputtering  yield.  E^ly 
experiments  on  single  crystals,  for  example,  suggested  that  the  yield  of  particles  is 
reduced  when  the  angle  of  incidence  occurs  along  a  channeling  direction  and  that 
it  is  increased  when  the  particle  is  incident  upon  a  random  direction  (Magnuson  tc 
Carlston,  1963;  Robinson  k  Southern,  1967). 

There  is  another  issue  relevant  to  this  discussion  which  is  embedded  in  the  dy- 
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Figure  13.  Parameter  <tcliiutioiu  for  the  ahaidow-coiie-enhaaced  deaorpti<m  experiment.  The  ion 
beam  i«  iitcidwt  at  tfi,  the  deaorbedpaiticlea  are  detected  at  and^  s  fi,  —$4.  The  thadoer-cane 
ia  described  hgr  a  radius  r  at  a  distance  I  behind  the  target  atom.  The  d  and  h  values  describe 
the  surface  bond  lengths. 


natnics  of  the  atomic  motion.  That  is,  specific  impact  points  on  the  surface  exhibit 
very  high  yields  while  other  points  are  rather  inactive  for  sputtering.  This  obser¬ 
vation  led  to  the  concept  of  atoms  sputtered  per  incident  ion  or  ASI  distributions 
(Harrison,  1988).  For  1  keV  Ar*^  ion  on  Cu,  for  example,  the  Cu  yield  varies  from 
0  to  over  15,  depending  on  the  impact  point.  It  would  be  interesting  to  be  able 
to  find  these  high  action  points  in  the  laboratory,  a  dim  prospect  considering  the 
aiming  accuracy  required. 

It  is  possible  to  take  advantage  of  the  above  information  to  utilise  the  angle 
of  incidence  of  the  ion  beam  to  glean  even  more  detailed  information  about  sur¬ 
face  structure.  The  collision  of  the  primary  ion  with  a  surface  atom  creates  a 
shadow-cone  which  focuses  the  incident  flux  to  specific  surface  coordinates  (Van 
der  Veen,  1985).  Shadow-cones  have  been  discussed  extensively  since  their  dis¬ 
covery  by  Lindhard  (1965).  As  seen  in  Figure  12  for  1  keV  Ar***  on  a  Ni  atom, 
this  flux  is  sharply  peaked  at  the  shadow-cone  edge  and  exhibits  a  width  of  less 
than  0.02  A.  When  the  tail  of  this  shadow-cone  strikes  a  surface  or  near-surface 
atom,  the  desorption  yield  should  dramatically  increase.  In  effect,  we  can  ‘rum’ 
the  incident  beam  to  specific  high  action  points  in  the  target.  With  knowledge  of 
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Flgur«14.  The  ilMdow-coMAnliaiiccd  SIMS  spectra  of  the  dcaorbedAf'^ -ion  jrieldM  A  ftinctioB  of 
the  angle  of  incidence.  The  particles  are  odlected  along  the  (a)  <100>,  (b)  <110>  and  (c)  <311> 
asinwith  during  3-keV  Ar^’-ion  bombardment  at  a  current  of  ~S  nA.  Only  particles  desorbed  in 
l^ane  with  0  *  35*  are  collected.  (Fbom  Chang  L  Winograd,  1989). 


the  shadow-cone  shape,  it  is  then  feasible  to  determine  surface  geometries  using 
simple  triangulation.  The  strategy  is  reminiscent  of  impact  collision  ion  scattering 
spectrometry  (ICISS)  where  it  is  possible  to  detect  the  crystal  orientation  which 
gives  rise  to  incident  ions  that  are  backscattered  by  180*  (Aono  et  al.,  1981;  Aono 
et  al.,  1983;  Yarmoffet  al.,  1986). 

The  geometry  of  the  shadow-cone  enhanced  desorption  experiment  is  illustrated 
in  Figure  13.  In  principle,  the  detector  should  be  configured  to  collect  all  desorbing 
ions,  regardless  of  their  kinetic  energy  or  angle.  Most  experiments  to  date,  however, 
have  been  performed  with  a  fixed  angle  0  and  with  a  detector  set  to  collect  only  the 
high-kinetic-energy  particles.  With  this  configuration,  the  shadow-cone  enhanced 
desorption  concepts  may  be  more  carefully  evaluated  (Chang  et  al.,  1987).  More¬ 
over,  it  has  been  shown  (Chang  ti  Winograd,  1989)  that  for  between  ±45*, 
ion  neutralization  effects  are  not  of  sufficient  magnitude  to  shift  the  angles  of  the 
intensity  maxima,  although  the  relative  position  of  the  peaks  may  be  altered  some¬ 
what.  This  is  a  convenient  geometry,  then,  for  indirectly  examining  the  behavior 
of  neutral  atom  ejection. 

Representative  results  (Chang  &  Winograd,  1989)  for  Ag{110)  are  shown  in 
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Figure  IS.  Ag-Ct  bond-length  dtangi  (O) ^nd  Auger  Cl  L3M3,3M}^  kinetic  energy  (•)  m  • 
function  of  CI3  expoeure  «t  300  K.  The  rfejAnted  value  refert  to  the  distance  between  the  center  of 
silver  atom  1  and  the  Cl  atom.  Thep(3xl)  LEED  pattern  was  observed  in  the  exposure  region  as 
shown.  The  bond  lengths  from  the  SEXAPS  experiments  are  associated  with  the  LEED  pattern; 
the  same  value  was  also  obtained  at  a  coverage  beyond  4  L  and  associated  with  the  c(4x2)  LEED 
pattern.  The  thadow-ctme  induced  desorption  mechanism  used  to  calculate  the  bond  length  is 
shown  in  the  inset.  The  interplanar  spacing  refers  to  the  relaxed  dean-surface  Ag{110}  value. 
The  shadow-cme  sluq>e  was  calculated  using  a  Thomas-Femu-Moliere  potential  with  a  scaling 
factor  of  0.86.  (FVom  Winograd  It  Chang,  1989). 


Figure  14.  There  are  two  major  features  in  these  distributions  which  correspond 
to  intersection  of  the  shadow-cone  with  a  second-layer  atom  {0  'v  35**)  and  with  a 
first-layer  atom  {0  ~  70”).  The  fact  that  these  peak  maxima  actually  correspond 
to  the  intersection  of  a  shadow-cone  edge  with  a  substrate  atom  has  been  tested  by 
using  a  full  three-dimensional  computer  simulation.  The  position  of  this  peak  has 
been  shown  to  be  unaffected  by  the  image  interaction  or  other  forces  related  to  the 
secondary  ionization  process  since  the  peak  position  is  found  to  be  independent  of 
0  for  9i  values  between  0”  and  45”.  From  the  position  of  the  peaks,  it  is  possible 
to  calculate  d  using  a  Moliere  potential  and  to  compare  this  calculated  value  to 
the  known  interatomic  spacing  for  Ag.  The  results  compare  to  within  ±  0.5”, 
yielding  an  uncertainty  of  about  ±0.06  A  in  the  bond  length.  Similar  procedures 
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for  deterinining  h  suggest  that  the  spacing  between  the  first  and  second  layer 
is  relaxed  by  (7.8  ±  2.5)%  and  the  spacing  between  the  first  and  third  layer  is 
relaxed  by  (4.1  ±  2.1)%  relative  to  the  bulk  spacing.  These  values  agree  within 
the  same  error  limits  with  those  bond  lengths  found  by  Rutherford  backscattering 
measurements  (Kuk  tc  Feldman,  1984). 

Adsorbate  atoms  will  also  create  shadow  cones  that  can  intersect  nearby  sub¬ 
strate  atoms,  opening  the  possibility  of  measuring  the  bond  length  of  chemisorbed 
species.  An  excellent  case  to  test  this  idea  is  the  chemisorption  of  Cl  on  Ag{110). 
Channeling  and  blocking  experiments  suggest  that  the  Cl  atom  binding  site  is  in¬ 
variant  over  a  wide  coverage  range,  allowing  the  Ag-Cl  bond  length  to  be  examined 
under  a  variety  of  experimental  conditions  (Chang  et  al.,  1987 ;  Winograd  &  Chang, 
1989).  The  results  for  the  shadow-cone  experiments  are  shown  in  Figure  15.  Of 
particular  note  is  that  at  extremely  low  Cl  coverages,  the  observed  bond  length  is 
extended  over  that  observed  in  the  high  coverage  limit  by  nearly  0.4  A.  This  change 
in  length,  accompanied  by  a  change  in  the  shape  of  the  Auger  electron  emission 
spectra,  has  been  explained  as  being  due  to  a  shift  from  highly  ionic  bonding  at 
low  coverage  to  more  covalent  bonding  at  high  coverage.  If  the  Ag-Cl  interaction 
is  highly  ionic,  the  bond  should  have  a  significant  dipole  moment.  Presumably 
dipole-dipole  interactions  between  nearby  Ci~  atoms  force  charge  back  into  the  Ag 
substrate.  The  absolute  value  of  the  bond  length  at  high  exposures  agrees  quite 
well  with  surface- EX  AFS  experiments  obtained  from  the  p(2x  1)  ordered  overlayer 
(Holmes  et  al.,  1987)  as  seen  in  Figure  15.  The  results  are  also  consistent  with  the 
channeling  and  blocking  experiments  discussed  previously,  which  suggest  that  the 
Cl  adsorbate  could  more  closely  approach  the  surface  at  higher  coverages  (Moon 
et  al.,  1986). 

These  types  of  surface  structure  determinations  are  readily  extended  to  more 
complex  systems.  Of  special  interest  are  compound  semiconductor  surfaces,  not 
only  because  of  their  technological  relevance,  but  also  because  they  display  a  rich 
diversitvof  surfacfcatructures  which  need  to  be  chuacterized.  A  good  example 
is  tiaAs{001}M2xx4)^hich  was  examined  earlier  in  this  review  with  regard  to 
its  behavior  unaSrUombardment  at  normal  incidence.  The  most  straightforward 
objective  of  the  shadow-cone  experiments  would  be  to  determine  the  length  of 
the  As- As  bond  in  the  surface  layer.  This  number  is  difficult  to  obtain  by  other 
methods,  but  is  important  for  input  into  calculations  aimed  toward  predicting  the 
surface  electronic  structure  (Chadi,  1987;  Larsen  &  Chadi,  1988). 

The  results  of  this  experiment  (Xu  et  al.,  1992)  are  shown  in  Figure  16.  The 
crystal  direction,  <011>,  corresponds  to  shooting  the  incident  beam  directly  along 
the  As2  bond,  in  an  azimuthal  direction  parallel  to  the  missing  row  of  As  atoms. 
The  spectrum  consists  of  four  major  peaks  at  polar  angles  of  0i  =  70.1*,  63.0®, 
44.5*  and  25.8*.  There  are  many  other  distinguishable  features  in  the  spectrum 
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Figure  16.  The  relative  inteneity  of  30-eV  Gt.'*'  iona  descurbed  firom  the  GaAa{001}(2x4)  euiface 
by  3-keV  Ar'*‘-ion  bombardment,  plotted  ae  a  function  of  the  ion-beam  polar  angle,  rdative  to 
the  lurface  normal.  The  ion  beam  is  parallel  to  the  <011>  crystal  direction.  (FVom  Xu  et  al., 
1992). 


which  have  been  confirmed  to  be  real  structure,  but  their  interpretation  awaits 
future  geneiations  of  research  students  (Burnham  et  al.,  private  communication). 
The  main  peaks  have  been  assigned  on  the  basis  of  parallel  molecular  dynamics 
calculations  for  Si  (Xu  et  al.,  1992).  A  summary  of  those  results  is  given  in  Table 
I.  Notice  that  the  As2  distance  is  determined  from  the  small  shoulder  at  63*  to 
be  2.73±  0.10  A.  The  procedure  for  calculating  the  bond-length  is  identical  to 
that  described  for  the  Ag{110}  surface.  This  value  of  2.73  ±  0.10  A  is  within 
experimental  error  of  a  grazing-incidence  x-ray  diffraction  value  from  GaA8{001} 
c(4x4)  of  2.59  ±A  (Sauvage-Simkin  et  al.,  1989).  These  numbers  suggest  that  the 
dimer  bond  distance  is  much  closer  to  the  bulk  As  bond  distance  with  threefold 
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Figure  17.  Rmtio  of  intensitica  (dN*  ldv/{dNldxj)  ve  1/vx  for  different  angle*  of  ejection.  The*e 
data  are  direct  ratio*  of  the  inteneitie*  given  in  Figure  1.  The  etraight  line*  have  been  fit  to  the 
high  velocity  portion  of  the  data  and  have  a  elope  of  A/a  (in  unit*  of  10*  cm/*)  a*  dieplayed  in 
eadi  frame  (FVom  Shapiro  k,  Fine,  1989). 


coordination  than  to  a  tetrahedral  covalent  As  bond  involving  sp^  hybridization 
(Kittel,  1986). 

The  fact  that  these  simple  angular  distributions  provide  such  microscopic  in¬ 
formation  about  surface  structure  is  really  remarkable  from  several  points  of  view. 
Certainly,  the  primary  beam  is  creating  lots  of  damage  to  the  surface,  and  yet  it 
seems  possible  to  obtain  accurate  surface  bond  lengths.  The  experimental  configu¬ 
ration  is  exceedingly  straightforward,  requiring  only  a  simple  ion  source,  quadrupole 
mass  filter  and  polar  angle  rotation  capabilities  for  the  sample  holder.  The  potential 
sensitivity  of  this  approach  to  low  concentrations  of  overlayers  is  indeed  unprece¬ 
dented.  The  three  examples  presented  in  this  review  are  very  promising  ones,  but 
much  research  remains  to  be  accomplished  in  this  area.  It  will  be  interesting  to 
examine  how  the  distributions  change  with  shadow-cone  radius  and  whether  it  will 
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be  possible  to  oihsnce  certain  desorption  mechanisms  by  changing  particle  mass 
and  energy.  Perhaps  the  most  important  aspect  of  this  experimental  configuration, 
however,  is  that  it  allows  an  accurate  description  of  the  sputtering  event  to  be  per¬ 
formed  without  using  molecular  dynamics  calculations.  The  angular  distributions 
can  be  largely  explained  using  only  the  details  of  the  first  encounter  of  the  primary 
ion  with  the  crystal  surface. 


6  Spectroscopic  Studies  of  Single  Crystal  Sput¬ 
tering 

The  laser  techniques  described  in  section  III  and  IV  are  inherently  state  selective  as 
illustrated  in  Figure  1.  For  example,  the  ground  state  can  be  ionized  through 
the  state  with  two  312.4  nm  photons.  The  metastable  ^Fr/a  state  lies  0.2 
eV  above  the  ground  state  and  can  be  ionised  using  two  328.0-nm  photons.  Since 
is  a  forbidden  transition,  it  is  possible  to  make  detailed  measurements 
on  both  states  and  to  compare  the  results  directly.  The  comparison  provides  an 
unprecedented  level  of  detail  regarding  the  mechanistic  aspects  of  energy  excitation 
and  quenching  at  surfaces. 

A  direct  measurement  of  the  complete  EARN  distribution  for  both  the  ^Fg/a 
ground  state  and  the  ^Fr/a  excited  state  for  Rh{100}  has  recently  been  completed 
(Winograd  et  al.,  1992;  Bernardo  et  al.,  1992),  and  both  maps  are  shown  in  Figure 
8.  Note  that  the  intensity  of  the  desorbed  atoms  in  their  ground  state  is  very  weak 
along  the  close-packed  direction  =  45*  as  defined  in  Figure  3)  due  to  blocking,  as 
expected  from  the  classical  calculations.  The  open  direction  (^  =  0*)  largely  arises 
from  second-layer  atoms  being  focused  upward  by  first-layer  atoms.  The  EARN 
map  of  the  ^F?/]  level  is  significantly  different  from  the  ground-state  distribution. 
There  is  considerably  more  relative  intensity  observed  at  normal  ejection  {B  = 
0*).  At  low  energies,  only  a  shoulder  exists  along  the  ^  =  45*  azimuth  (<110> 
direction).  And  finally,  the  falloff  with  energy  is  much  slower  for  the  excited  state, 
most  likely  that  the  population  of  the  *Fi/2  state  is  a  direct  measure  of  the 
/^excite^robability,  uncomplicated  by  so-called  cascading  effects.  This  conclusion 
/  T^deriTCd  from  two  important  facts.  First,  the  intensity  of  the  ground-state  signal 
is  approximately  20  times  the  signal  observed  for  the  *F7/2  state.  Second,  the 
C'f-Cihurii^  intensity  of  the  next  higher-lying  excited  state  (^Fs/2  with  excitation  energy  of 

0.3  eV  above  the  ground  state)  is  at  least  2  orders  of  magnitude  smaller  than  the 
^F7/2  state.  Earlier  studies  by  Young  and  co-workers  (1984)  using  laser-induced 
fluorescence  have  shown  that  the  population  density  for  Fe  decreases  exponentially 
with  the  magnitude  of  the  excited  energy,  roughly  in  accord  with  the  observations 
for  Rh.  Similar  results  have  been  obtained  using  MPRI  (Kimock  et  al.,  1984). 
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Hmce,  time  ue  en  inaignificant  number  of  excited  Rh  atoma  in  aufSciently  high- 
lying  excited  atatea  to  contribute  to  the  ^/V/j  level  by  caacading  from  higher  levela. 
Once  theae  atoma  reach  the  gaa  phaae,  then,  it  is  unlikely  that  the  population 
denaity  ia  affected  by  decay  from  other  excited  atatea. 

The  interpretation  of  theae  energy  and  angle-reaolved  excitation  probabilitiea 
are  only  now  being  aorted  out.  The  general  idea  ia  to  try  to  describe  the  dynam- 
ica  of  the  excitation  and  deexcitation  physica  using  known  quantum  mechanical 
phenomenon  and  then  to  examine  the  effects  of  this  dynamics  when  it  ia  coupled 
together  with  classical  molecular  dynamics  (Bernardo  et  al.,  1992).  Sevnal  im¬ 
portant  concepts  are  emerging  from  these  models.  First,  the  form  of  the  velocity 
dependence  is  strongly  reminiacent  of  Hagatrun^  (1954)  original  energy  transfer 
ideas.  A  model  which  ia  consistent  with  this  picture  is  currently  being  exploited 
by  several  groups  and  involves  collisional  excitation  followed  by  nonradiative  de¬ 
excitation  (Bernardo  et  al.,  1992;  Shapiro  it  Fine,  1989).  Both  the  excitation  and 
deexcitation  steps  obviously  depend  upon  the  nature  of  the  collision  cascade  as 
it  evolves  in  the  crystal  and  on  the  electronic  properties  of  the  material.  For  the 
data  shown  in  Figure  17,  note  that  the  slopes  of  the  change  in  the  logarithm  of 
the  excited  state  fraction  with  1/vj,  depends  on  both  polar  and  aaimuthal  angle. 
This  effect  is  not  predicted  without  including  the  influence  of  the  cascade  on  the 
entire  process  (Winograd  et  al.,  1992;  Bernardo  et  al.,  1992;  Winograd,  1992).  For 
example,  if  glancing  collisions  are  important  at  high  polar  angles,  the  quenching 
rate  will  be  appropriately  higher. 

Although  there  are  many  potential  intricacies  associated  with  these  experi¬ 
ments,  it  is  of  interest  here  to  conclude  with  one  final  observation.  There  is  a 
region  of  the  data  shown  in  Figure  17  where  the  excited  fraction  is  found  to  be 
independent  of  velocity.  Although  a  number  of  factors  contribute  to  this  effect,  the 
molecular  dynamics  calculations  show  that  collisions  of  atoms  over  the  surface  are 
most  important  (Winograd  et  al.,  1992;  Bernardo  et  al.,  1992;  Winograd,  1992). 
With  this  mechanism,  two  independently  ejected  atoms  undergo  a  hard  collision 
over  the  surface  and  become  reexcited.  Because  they  are  out  of  the  electronic 
coupling  range  of  the  solid,  the  excitations  are  not  quenched. 


7  Conclusion  and  Prospects 

The  revolution  in  modern  surface  science  has  opened  new  opportunities  for  learning 
more  detaib  about  single  crystal  sputtering  than  ever  before.  Moreover,  new  de¬ 
tection  schemes  have  made  it  possible  to  essentially  examine  the  collision  dynamics 
on  an  event  by  event  basis.  The  new  data  are  now  in  quantitative  agreement  with 
results  from  sophisticated  molecular  dynamics  cdculations,  at  least  for  a  limited 
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number  of  model  systems.  This  agreement  allows  one  to  extract  the  important  en¬ 
ergy  dissipation  mechanisms  Irom  the  computer  printouts  and  to  have  confidence 
in  their  validity. 

Yet,  there  is  much  to  be  done.  The  computer  simulations  are  complex  and  the 
experhnentalists  need  simpler  formulas  by  which  to  interpret  their  data.  These 
formulae,  however,  should  provide  a  quantitative  prediction  of  the  sputtering  yield 
as  a  function  ejection  energy  and  angle.  Perhaps  the  shadow-cone  ^;>ptoach, 
where  the  mathematics  is  easy,  is  a  step  in  that  direction.  Spectroscopic  studies 
of  single  crystal  sputtering  are  only  beginning  and  offer  hope  of  disentangling  the 
factors  that  go  into  inelastic  energy  loss.  For  example,  the  irtiiflrn  diirnwH  hrrr 
invdved  only  electronic  states  that  are  part  of  the  ground  stat^monofijj^igher 
lying  states  will  undoubtedly  yield  new  surprises.  And  finally,  there  nei^  to  be 
mote  wwk  on  molecular  crystals,  allojrs  or  on  insulating  crystals  in  order  to  continue 
to  elucidate  the  mechanistic  differences  associated  with  the  sputtering  of  these 
disparate  materials.  Spectroscopic  studies  of  neutrd  clusters  emitted  from  single¬ 
crystals  under  low-dose  conditions  would  provide  interesting  information  about 
the  vibrational  and  rotational  excitations  that  occur  in  these  species.  There  will 
certainly  be  many  interesting  new  effects  discovered  with  these  substances.  Their 
role  is  already  becoming  increasingly  important  in  modern  materials  science. 
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Crystal  Direction 

Peak 

Dbulk 

D 

<011> 

Peak  1 

70.1 

4.00 

4.G0  ±  0.10 

Dm 

Peak  2 

63.0 

4.00 

2.73  ±  0.10 

Du 

Peak  3 

25.8 

1.41 

1.40  ±  0.10 

Djs 

Table  I.  Experimental  and  calculated  results.  The  distances  are  in  unit  of 
Angstroms  and  the  angles  are  expressed  in  degrees.  Du  u  the  first-layer  As-As 
bond  distance  in  a  <01 1>  direction;  Dn  the  second-layer  Ga-Ga  distance  along  a 
<011>  direction. 


